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Abstract. In this talk, I address the comparison between results from lattice QCD 
computations and Chiral Perturbation Theory (ChPT). I briefly discuss how ChPT 
can be adapted to the much-used quenched approximation and what it tells us about 
the special role of the rf in the quenched theory. I then review lattice results for some 
quantities (the pion mass, pion scattering lengths and the K + — > -k + tx matrix element) 
and what quenched ChPT has to say about them. 



1 Introduction 

Chiral Perturbation Theory (ChPT) (Weinberg (1979), Gasser and Leutwyler 
(1985), Gasser (1997)) gives us information about the functional dependence 
of quantities associated with low-energy Goldstone Boson (GB) physics on the 
light-quark masses. Examples of such quantities are the GB masses, decay con- 
stants and scattering amplitudes. At any given order, these relations between 
physical quantities and quark masses involve a finite number of constants (the 
"low-energy constants (LECs)," which cannot be determined from ChPT alone. 
Therefore, ChPT is predictive (at any given order) if we consider a number of 
physical quantities larger than the number of LECs needed at that order. These 
LECs can in principle be determined either by comparison with experimental 
data, or by a theoretical calculation from the underlying theory, the Standard 
Model. 

The strong-interaction part of such calculations is nonperturbative, and this 
is where Lattice QCD (LQCD) comes in. In LQCD, physical quantities (or re- 
lated quantities, such as weak matrix elements) are computed from first princi- 
ples as a function of the quark masses. By fitting the results with the relations 
predicted by ChPT, one can then, in principle, determine the LECs. This is very 
similar to determining these constants from experimental data, with the added 
advantage that in LQCD one can vary the quark masses. 

Of course, in general, LQCD results will need to be of a high precision in or- 
der to extract the 0(p 4 ) LECs, because, in general, they show up in the one-loop 
corrections to the tree-level predictions from ChPT. This means good control 
over both statistical and systematic errors in the lattice computations. For in- 
stance, we need the volume 1? in lattice units to be large in order to use a 
small lattice spacing a, while keeping the physical volume large enough to fit the 
hadronic system of interest. It is important to keep in mind that lattice results 
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need to be extrapolated to the continuum limit before they can be compared with 
ChPT. Also, for ChPT to be valid, GB masses need to be small compared to 
the chiral symmetry breaking scale. Again, this leads to the requirement of large 
enough volumes. (Finite volume effects can be studied within ChPT (Gasser and 
Leutwyler (1987)); we will see some examples in the following. However, LQCD 
practitioners are not only interested in the comparison with ChPT!) The tempo- 
ral extent of the lattice needs to be large enough so that the lowest state in any 
channel can be reliably projected out by taking the large-time limit of euclidean 
correlation functions. 

This has led to the widespread use of the Quenched Approximation (QA) 
(Parisi et al. (1981), Weingarten (1982)), in which the quark determinant is 
omitted from the LQCD path-integral. This is equivalent to omitting all contri- 
butions to correlation functions that involve sea-quark loops (valence-quark loops 
resulting from contractions of quarks in composite operators of course remain 
present, as they have nothing to do with the determinant). The reason is that, if 
the fcrmion determinant is included, the necessary computer time increases by 
orders of magnitude if we keep the physical parameters (lattice spacing, volume, 
quark masses) the same. 

Quenched QCD is a different theory from QCD (as we will see, it's not even 
a healthy theory!), and therefore the predictions of ChPT do not apply to lattice 
results obtained in the QA. Fortunately, it turns out that a quenched version of 
ChPT (QChPT) can be developed systematically, and it is in this framework that 
we can compare quenched LQCD with ChPT. Unfortunately, this also implies 
that the LECs predicted by quenched QCD are not necessarily equal to those 
of full QCD, and as such, knowledge of them is of somewhat limited value. In 
this talk, I will describe how QChPT works, and discuss some examples of the 
comparison of lattice results with ChPT. I should note right away that, as we 
will see, currently lattice data are not precise enough yet to unambiguously see 
ChPT one-loop effects. 

Before we get into this, let me end this introduction with a few remarks. First, 
recently, more lattice results with "dynamical fermions" (i.e. including sea-quark 
loops) are becoming available. However, since the overhead in computing the 
fcrmion determinant is so large, these results are often for one or two values of the 
sea-quark mass, while many values of the valence quark masses are considered. 
The methods described in this talk can easily be adapted to these "partially 
quenched" theories with sea quarks with a mass that differs from that of the 
valence quarks (Bernard and Golterman (1994), Sharpe (1997b)). Second, the 
study of QChPT sheds much light on the nature of the QA, and, as such, has 
been very helpful for LQCD. 

2 Quenched ChPT, the rj', and the Pion Mass 

Euclidean quenched QCD can be defined by taking the ordinary QCD lagrangian, 
and adding a new set of quarks {^} to it which carry one- by-one exactly the 
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same quantum numbers as the normal quarks {qi}, but which have opposite 
statistics (Morel (1987)): 

£ = Qi(P + m%)<li + q l (P + m i )q l , i = u,d,s. (1) 

We will refer to these new wrong-statistics quarks as "ghost-quarks." The reason 
for their introduction is that the contribution coming from integrating over these 
ghost-quarks exactly cancels the determinant which comes from the integration 
over the normal quarks. It is easy to see this diagrammatically: while a minus 
sign is needed for every occurence of a quark loop, this is not the case for a 
ghost-quark loop, because of their bosonic statistics. So, for every diagram with 
a quark loop, a similar diagram contributes with that one loop replaced by a 
ghost-quark loop, therefore leading to the cancellation of this diagram. 

For mi = 0, the lagrangian (1) is invariant under a six- flavor chiral symmetry 
group (Bernard and Golterman (1992)). A typical transformation looks like 

(q' i \ = (A ij Bij \ f Qj\ 

Uv yea aJUJ 

where A, B, C and D are 3x3 blocks, with B and C containing anticommuting 
numbers, since they transform bosons into fermions and vice versa. This group 
is known as a graded version of U(6), and is denoted by U(3|3). For rrii — 0, the 
symmetry of (1) is therefore U(3|3) L xU(3|3) fl . 

We will assume that gluons are responsible for forming mesonic bound states. 
In this case, since gluons couple equally to quarks and ghost-quarks, that means 
that not only the usual qq GBs (i.e. it, K, w and rj') will occur, but also qq "ghost- 
mesons," and qq, qq fermionic "hybrid" mesons. (For the role of rj', see below.) 
QChPT can then be developed systematically in a way completely analogous to 
the usual three-flavor case (Bernard and Golterman (1992), for other early work 
on QChPT, see Sharpe (1990,1992)), but now for the group U(3|3) L xU(3|3) fl . 
The quenched effective lagrangian describes the low-energy physics of all Gold- 
stone mesons, ghosts and fermionic hybrids included. 

Lack of space prevents me from giving more details on the technicalities of 
QChPT (see e.g. Golterman (1994)). However, before we continue to consider 
any examples, we need to address the very special role of the rf . In unqucnchcd 
QCD, the r/' is not a GB because of the U(1)a anomaly. In quenched QCD, we 
would expect that a similar role is played by some linear combination of 77' and 
fj'. We can understand which linear combination by realizing that ghost-quark 
loops do not contribute a minus sign, so that we get a nonzero triangle anomaly 
between the U(l) axial current and two gluons by subtracting the triangle graphs 
with a ghost-quark running around the loop from those with normal quarks. 
Hence in the quenched theory, the "anomalous" meson is the "super-?/" , (rj' — 
r/')/\/2. The super-77' transforms as a singlet under the nonanomalous part of 
the graded chiral symmetry group, and therefore arbitrary functions of this field 
(and its derivatives) appear in the effective lagrangian (Gasser and Leutwyler 
(1985), Bernard and Golterman (1992)). 



(2) 
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Let us consider the lowest order, 0(p 2 ) part of the effective lagrangian 
quadratic in the rj' and ff fields: 



= M) 2 -MY + ^mUv' 2 V' 2 ) (3) 



+ln 2 (v'-v') 2 + l*(Wv'-v')) 2 , 



where the terms on the second line are the quadratic, 0(p°) and 0(p 2 ) parts of 
the arbitrary functions of the super-?/ field. All minus signs originate from the 
graded nature of the chiral symmetry group (in building invariants for such a 
group, the supertrace plays the role of the trace in the ordinary case, see De- 
Witt (1984)). They will have important consequences. For simplicity we assumed 
degenerate quark masses. 

If we omit the rj' field from (3), we obtain the quadratic part of C 7 y for the 
unquenched theory, and the parameter \j? would essentially be the singlet part 
of the rf mass (one can show that ji 2 does not vanish because of the anomaly!). 
So, let us see what happens in the quenched case. We will treat the first line in 
(3) as the part defining rf and ff propagators S®, (p) and S~, (p) : 

S°(p) = -^,(p) = _-L_. (4) 

The second line defines the two-point vertices — {ji 2 + ap 2 ) for r)'-r]', ff-ff and 
for rf-ff mixing. We can find the complete rf two-point function by summing 
all diagrams containing an arbitary number of these two-point vertices on an rf 
line (taking into account combinatoric factors): 

S n ,(p) = S v ,(p) (5) 
-S v ,(p)(» 2 +ap 2 )S°,(p) 

+S»(,i 2 + ap 2 )(S° n ,(p) + S°,(p))(p 2 + ap 2 )S°,(p) 



Using (4), we see that the third line in (5) vanishes (and so do all higher order 
contributions), leading to 

1 ii 2 ~\~ otp 2 

Sv ' {P) = WT^ 2 ~ {p 2 +ml) 2 ■ (6) 

It is actually easy to convince oneself, that this cancellation is nothing else than 
the cancellation of sea-quark loops with their ghost counterparts. The vertex 
— (fi 2 + ap 2 ) represents the annihilation of the valence quark antiquark pair in the 
r}' (which is a flavor singlet, making this annihilation possible), and the creation 
of either a quark antiquark, or a ghost-quark ghost- antiquark pair. Tracing the 
quark flow for each of the terms in (5), one finds that quark and ghost-quark 
loops occur for all terms with more than one vertex insertion. 
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From (6), it is clear that the rj' is "sick" in the QA because there is no particle 
interpretation for the double pole term. Moreover, it tells us that we have to keep 
the r\' in the QA, because it has poles which are degenerate with those of the 
other Goldstone mesons. These poles can lead to chiral logarithms that have 
no counterpart in the unquenched theory. In the case of nondegenerate quark 
masses, the n° and 77 will inherit the problems of the rj' through mixing. 

As an example, consider the one-loop expression for the pion mass in terms 
of the quark mass, calculated in QChPT: 



where, for simplicity, we set a = (of course in comparing with real data, one 
has to take a into account). A, B and C are parameters related to the LECs, 
which are not predicted by ChPT. We see that the chiral log in (7) is entirely due 
to the double pole term in the rj propagator (because it is proportional to /i 2 ). 
In unquenched ChPT, the chiral log would be multiplied by an extra factor m q , 
and be more suppressed for small quark mass. In fact, in the quenched case, the 
quantity mj/m q diverges in the chiral limit, m q — > 0! This is the first example 
of the highly infrared divergent behavior of quenched QCD. (I believe that this 
divergence is a sickness of quenched QCD, and not just of QChPT (Bernard and 
Golterman (1993)).) 

In Fig. 1, I show lattice results for this quantity (in lattice units), from a 
compilation of staggered fermion spectrum results by Gottlieb (1997) (see this 
review for the original refs. on all these data). The mass ratio on the horizontal 
axis is basically the quark mass. (3 — 5.7 corresponds to a lattice spacing a ~ 
0.2 fm, while [3 = 6.5 corresponds to a w 0.055 fm. The squares and fancy 
crosses correspond to the data with the largest spatial volumes. The overlapping 
squares have L 2.6 fm and L « 3.5 fm, while the highest square (for lowest 
quark mass) has L « 1.8 fm. The fancy crosses have L m 2.6 fm. The fact that 
the squares for L w 2.6 fm and L s=a 3.5 fm overlap indicates that L w 2.6 fm 
is large enough to not have finite volume effects larger than the size of the 
error bars. Points for different (3 do not fall onto one curve because the quantity 
on the vertical axis depends on a. While the data do show an upward trend 
with decreasing quark mass, as predicted by (7), we cannot conclude that this 
behavior is unambiguously seen in these data. We know that scaling violations 
occur for the data at the two lower (3 values, so that these cannot be reliably 
compared to the continuum prediction of QChPT. At larger (3, the curves look 
flatter, and more data for small quark masses at these higher values of (3 will be 
needed to fit (7). Qualitatively, we do see a departure from the ChPT tree- level 
prediction m 2 oc m q , which would correspond to a horizontal line in Fig. 1. For 
a recent, much more detailed discussion of attempts to fit (7) to lattice data, 
including a list of subtleties and pitfalls, see Sharpc (1997a). We note here that it 
is probably better to determine 5 from a certain ratio of decay constants (Bernard 
and Golterman (1993), Bhattacharya and Gupta (1996), Sharpe (1997a)). 




Arriq (1 + S log Bm q + Cm q ) 
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Fig. 1. Quenched staggered fermion lattice results for m^-/m q , in lattice units, versus 
the quark mass in dimensionless units (see text). From the Lattice 96 spectrum review 
(Gottlieb (1997)). 



The logarithm in (7) appears with the coefficient 6, instead of the usual 
/ (4n f n ) 2 , and therefore is of the same order as the tree- level term in chiral 
power counting. At more than one loop, there are contributions of the same order 
in m^/(47r/ 7r ) 2 , but higher order in 5. This means that, in order for QChPT to 
be systematic, S has to be treated as a small parameter. In unquenched QCD, we 
can estimate 5 from the rf mass, yielding S ~ 0.18. Of course, since the quenched 
theory is different, S could have a different value in this case, and it is important 
to determine its value numerically. It is believed that the quenched value docs 
not differ much from the unquenched value, but errors are not yet sufficiently 
under control to quote a number (Sharpe (1997a)). We also note here that both 
^ 2 and a are of order 1/N in the large- AT expansion (Veneziano (1979), Witten 
(1979)). 

There is an extensive body of work on QChPT, and I do not have space here 
to properly review even a substantial fraction of it. There has been work on the 
inclusion of other hadrons (baryons: Labrenz and Sharpe (1996), vector and ten- 
sor mesons: Booth et al. (1997), Chow and Rey (1997), heavy baryons: Chiladze 
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Fig. 2. Relative error from quenching in /s, /s s , Bb and _Bb s at one-loop in QChPT. 
For explanation and assumptions, see text. From Sharpe and Zhang (1996). 



(1997)), weak matrix elements (Bk- Sharpe (1992), fs and Bb'- Booth et al. 

(1995) , Sharpe and Zhang (1996), K + —* 7r + 7r° and Bk'- Golterman and Leung 
(1997), baryon axial charge: Kim and Kim (1996)), pion scattering (Bernard and 
Golterman (1996)) and formalism (Colangclo and Pallantc (1997)). For work on 
fitting lattice data to results from QChPT, I refer to the Proceedings of the Lat- 
tice 96 and Lattice 97 conferences, in particular the reviews by Gottlieb (1997), 
Sharpe (1997a) and Yoshie (1997), and references therein. 

Here, let me show a sample result of a case where the GBs are coupled to 
other hadrons, in this case the B meson. Fig. 2, taken from Sharpe and Zhang 

(1996) , shows the relative error from quenching for fs, /b s , Bb and Bb s , as 
a function of the many parameters that are not determined by ChPT, from 
a one-loop ChPT calculation. Note that more parameters typically show up 
than in the case of pure GB physics, because of the various coupling constants 
between the GBs and the heavy sector, g (go) is the £?-pion coupling in the 
unquenched (quenched) theory, g' is the B-rj' coupling (which is suppressed in 
l/N). The figure assumes gQ — g, and values as shown (the value for g is from 
an estimate in the unquenched theory). All 0(p A ) LECs have arbitrarily been set 
to zero, while the cutoff is chosen at A = 1 GeV. Fig. 2 shows that, under these 
assumptions, quenching errors in these quantities can be large. As such, it gives 
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us useful information, but it should be stressed that this is not a calculation of 
quenching errors, because of the many possibly unjustified choices made for all 
the parameters. This situation is typical in the case that we consider the effective 
theory for GBs coupled to other hadrons. 

3 Pion Scattering 

LQCD computations give us access to euclidean correlation functions only, and, 
lacking the possibility to analytically continue them to Minkowski space, it ap- 
pears impossible to extract information on scattering lengths from the lattice. 
However, we can make use of the fact that lattice computations are necessarily 
done in a finite volume, and use this as a probe of the interactions between two 
pions. More precisely, one can compute the energy of a state with two pions 
at rest from the euclidean correlation function (for I = 2, a similar correlation 
function can be defined for 1 = 0) 

Ci =2 = (0|7r + (t)7r + (t)7r-(0)7r-(0)|0> (9) 
= ^e-^ t |H 7 r-(0)7r-(0)|0)| 2 

= Zc' 2m ^ (l-AEt + 0(t 2 )) + excited states, 

where 7r _ (0) creates a zero-momentum ir + at time 0, and 7r + (i) annihilates a 
zero- momentum 7r+ at time t. We inserted a complete set of states {|a)}, and 
assume that we can take t large enough to project out the lowest-energy state. 
This is the state with two pions at rest in finite volume, with = 2m w + AE, 
where AE is the energy shift due to finite volume effects (in infinite volume, two 
pions at rest do not interact). We expanded in AE t, because it is in this form 
that the energy shift will show up in a ChPT calculation. Also, it was this form 
which has been fitted to lattice data by various groups in attempts to determine 
AE (Guagnelli et al. (1990), Gupta et al. (1992), Kuramashi et al. (1993)). From 
(9) we see that AE t has to be small enough for the expansion on the last line 
to be valid (but of course t has to be large enough to project out the excited 
states). 

As mentioned above, AE is a measure of the interactions between the two 
pions, due to their confinement to a finite volume. A precise connection between 
AE and the pion scattering length was given by Liischer (1986): 

ae - -S? C 1 - 2 - 83ra97 f + 6 ' 375i 4) + ° G?) • < io » 

where ao is the infinite volume scattering length for the corresponding channel 
(7 = or 2). While the proof of (10) is more general, the r.h.s. can be recovered 
term by term at tree level, one loop, two loops, etc. in ChPT. It is instructive to 
look at this in some detail, because it will tell us how things change in the QA. 
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Consider the one-loop diagram in which the pions are created at time 0, 
scatter at time t 2 > 0, rescatter at time t\ > t 2l and then are annihilated at 
time t > t\. (There arc many other one- loop contributions of course, but it is 
this one that leads to the 1/L 4 term in (10).) Since the initial and final pions 
have zero spatial momentum, their propagators are of the form exp(— m n t) if 
they "travel" a time t; the expression for this diagram looks like 

A 2 jL 3 e - 2m „ t J_^ e -2(B(k)-m„)(t 1 -t 2 ) ; £( k ) = ^m 2 + k 2 . (11) 
k 

The second exponential takes into account that between the two scatterings, the 
two pions can have arbitrary momenta k and — k. The overall factor L 3 comes 
from integrating the diagram over the spatial volume, and A oc m 2 / (4n f^) 2 is 
the interaction strength. This expression has to be integrated over t\ and t 2 , 
where here we are interested in the contribution with < t 2 < t\ < t. This 
integration leads to 

* ^5^)-™- (12) 

(the k = term contributes to the 0(t 2 ) term in (9)). In infinite volume, we 
may replace X)k^o ^ v / d 3 k, and (after renormalization) we get a contribution 
~ A 2 /L 3 to AE. (The L 3 in (11) turns into a 1/L 3 after dividing by the factor L 6 
which normalizes the contribution L 6 exp(—2m 7T t) of two noninteracting pions.) 
This just gives us a one-loop correction to ao in the 1/L 3 term in (10). In a 
finite volume, the smallest momenta k ~ 1/L cause the sum to deviate from 
the integral most prominently. Expanding E(k) — <~ 1/L 2 (keeping track of 
the L dependence only), the sum (12) gives a contribution ~ t/L. This leads to 
an additional finite volume contribution to AE of order A 2 /L 4 ~ a^/L 4 , cf. the 
second term in (10). 

We can now also see how quenching, and in particular, the special role of 
the r]', will change this result. In the 1 = channel, the intermediate mesons 
in the diagram I just discussed can also be singlets, and the two-point vertex 
— (p 2 + ap 2 ) can appear on these internal lines. Let us consider the case of one 
such insertion at time tx- The expression for the diagram is similar to (11), 
but now we have to integrate also over t x ] - Considering the contribution where 
t 2 < tx < h, we pick up an extra factor t\ — 1 2 after doing the integral over tx, 
and from the integration over t\ and t 2 we now get 



1 \ - 



(13) 



L3^ (jB(k) _ m7r) 2 • 

Going through the same argument as before, this leads to an extra factor L 2 in 
the finite- volume energy shift, i.e. a contribution ~ S/L 2 to AE, quite unlike 
the prediction of Liischer's formula! Even worse, if we insert the r/' two-point 
vertex on both internal lines, we find another factor L 2 enhancement, and a one- 
loop contribution to AE of order S 2 /L° = S 2 . These "enhanced finite volume" 
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Fig. 3. Pion scattering lengths as a function of m n /m p . Errors are statistical only. 
The squares and diamonds are from Gupta et al. (1992), all other points are from 
Kuramashi et al. (1993), from which this figure was taken. 



corrections are another example of the serious infrared problems that affect the 
quenched theory, due to the double pole in (6). For a complete analysis (and real 
calculation), see Bernard and Golterman (1996). We just note here that also the 
rest of the structure of (10) is not reproduced in the quenched theory. For 1 = 2 
the problem is less serious, because no enhanced corrections occur (the internal 
lines cannot be singlets in this case). 

We conclude that QChPT does not only break down for small L, but also for 
large L, since the tree- level term in AE goes like 1/L 3 , while there are one- loop 
corrections which go like a lower power of 1/L. In fact, one has to worry whether 
QChPT applies for any choice of parameter values in this case! It is therefore 
very interesting to compare these results with quenched numerical computations 
of the pion scattering lengths. 

In Fig. 3, which I took from Kuramashi et al. (1993), lattice results for Oq -0 ' 2 
are shown, in units of tree-level ChPT predictions, as a function of m^/irip. Also 
results from Gupta et al. (1992) for I = 2 are included (squares and diamonds). 
The volume is 12 3 x 20 in lattice units, (5 = 5.7, corresponding to o k 0.2 fm. 
We see that, within errors (only statistical errors are shown), the lattice results 
agree with ChPT at tree level. In Bernard and Golterman (1996) we considered 
the size of the one-loop QChPT corrections for the parameters of these lattice 
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computations, and found that they are always smaller than 20% of the tree-level 
values. However, we cannot directly compare the lattice results with ChPT be- 
cause of the presence of other systematic errors, most notably scaling violations 
(a is too big) , and contamination of excited states (t in (9) may not have been 
large enough to reliably project out the lowest state, see Bernard and Golterman 
(1996)). 

4 K — ► tt+tt Decay 

LQCD practitioners have had a long-standing interest in computing nonleptonic 
K-dec&y matrix elements, because of such prominent experimental observations 
as the AI = 1/2 rule, and their importance in determining Standard Model 
parameters (the CKM-matrix). I will restrict myself here to K + — ► 7r + 7r°. Refer- 
ences Gavela et al. (1988), Bernard and Soni (1989) report on the status of past 
lattice computations of (7r + 7r° IO4 , while Ishizuka et al. (1997) reports on a 
very recent computation. Here O4 = (slI^^^lI^ul) + (sLj^u^iuLl^dL) — 
{~SLl^d L )(d L ^d L ) is the AI = 3/2, AS = 1 component of the SU(3) L 27-plet 
responsible for this decay. In ChPT, we have, with S the nonlinear GB-ficld, 

O4 = a^Sd^^Sd^S^)} + 0(p 4 ) operators, (14) 

where tl> projects out the appropriate component, a 27 is a free parameter, and 
the 0(p A ) operators were listed in Kambor et al. (1990). In Golterman and 
Leung (1997), (7T + 7r°|04|.ft' + ) was calculated to one loop in order to compare 
the real world with lattice results. Since the 0(p A ) coefficients are not or very 
poorly known, we set them equal to zero, and estimated the error from this by 
considering two values of the cutoff A = 770 MeV/1 GeV. The real- world value 
then is (for tree level, see Donoghue et al. (1982), for one loop with — 0, see 
Bijnens et al. (1984)) 

, + n,^ , I2ia 27 , , 2 N / 0.63,-4 = 1 GcV\ , , 

tt+tt \Oi \K + ) = ' (mi - ml) x 1 + n ' t ). (15 

x V 0.36, yl = 770 MeVy v ; 

On the lattice, one computes the matrix element from 

C(t 2 ,h) - (0|7r+(t 2 )7r°(t 2 ) 4 (h) K-(0)\0) (16) 

*2»*l>0 c -E 27 ,(t 2 -t 1 ) c -m K t 1 

(0|7T+ (0)7T°(0) Itt+tt ) (7T+7r°|C' 4 (0)|i^+) (K + \K~ (0) |0) 
X (ir+ir°\ir+ir°)(K+\K+) 

where all mesons are taken to be at rest. In actual lattice computations, degen- 
erate quark masses were used, so that ttlk = m w . The matrix element in (16) 
therefore is not the one we want: it is unphysical because of the choice of masses 
and external momenta. Furthermore, there are power-like finite volume effects, 
and the QA has been used in all computations. All these systematic effects can 



12 Maarten Golterman 



be estimated in one-loop ChPT (the unphysical choice of masses and momenta 
already at tree level in ChPT (Bernard and Soni (1989)). 

To one loop in ChPT the physical matrix element and the unphysical quenched 
lattice matrix element (after extrapolation to the continuum limit) are related 
by (Golterman and Leung (1997)) 



(7r+^ |O 4 (0)|X + ) 



with 



phys 



Y 



Y 



quenched 



2M1 



* (^^104(0)1^)™; 



0.089, 
-0.015, 



4=1 GeV 
4=770 MeV 



1 + 



Ml 



-3 log ^+F{M„L) 
17.827/(m OT L) + 12f 2 /KI) 3 . 



F(m„L) 



(17) 

(18) 
(19) 



Here parameters with a subscript q are those of the quenched theory, and M T = 
Mk, F v refer to the lattice quantities. At tree level, we would have Y = 1. 0(p 4 ) 
LECs have been set to zero. Some remarks are in order: 

• For degenerate quark masses, O4 does not couple to the 7/, which explains 
why (18) does not depend on S and a. For nondegenerate quark masses, such 
dependence would show up, presumably making the QA less interesting in 
that case. 

• The ratios a 27 /a% 7 and f q /f are not known, and, in the following discussion 
of lattice results, we will arbitrarily set them equal to one. 

• The 0(p 4 ) LECs are basically "absorbed" into the cutoff A. As before, we will 
take the cutoff equal to 770 MeV or 1 GeV independently in the quenched and 
unquenched theories, and take the spread as an indication of the systematic 
error associated with our lack of knowledge of these LECs. 

In Fig. 4 I show lattice data for the matrix element from Bernard and Soni 
(1989) (open symbols), where the tree-level correction factor {m 2 K — m^)/2M% 
in (17) was already taken into account. The errors on these points are statistical 
only, and we left out points at pion masses with M n > m p and/or at smaller vol- 
umes. The crosses show what we obtain if we "correct" the central values of each 
of these points with the factor Y (evaluated at the appropriate pion mass and 
volume). The errors on these points indicate the spread from choosing different 
combinations of A and A q . We see that at all points Y < 1, and that therefore 
the one-loop corrections reduce the discrepancy between lattice data and the 
experimental result. However, since one-loop effects are rather substantial, two- 
loop effects can probably not be neglected. Also, again, scaling violations and 
various other systematic effects are not taken into account. For a more complete 
discussion, see Golterman and Leung (1997). 

Very recently, the computation was done again by the JLQCD collaboration, 
with larger volumes, and larger (3 (Ishizuka et al. (1997)). The results are shown 
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Fig. 4. if + — ► 7T + 7r° decay amplitude as a function of pion mass. Open symbols: 
data from Bernard and Soni (1989) (squares: 16 3 x 25 (or x33), (3 = 5.7; oc- 
tagons: 24 3 x 40, [3 = 6); crosses: including the correction factor Y. The constant 
c = 2^/2/ (G F sin C cos 6» c ). 



in Fig. 5 (without the one-loop correction factor Y, along with the data of 
Bernard and Soni (1989)) and Fig. 6 (with the factor Y taken into account). 
If the unquenched cutoff A = A cont = 1 GeV is chosen, the values in Fig. 6 
would come out about 10% higher. Again, I should stress that, because of the 
uncertainties in the estimates for Y, the lack of knowledge of the ratios a 27 /ai 7 
and fq/f, and various systematic effects which cannot be estimated in ChPT, 
one can only conclude that one-loop ChPT reduces the discrepancy between 
lattice and experiment. 



5 Conclusion 

It is clear that one-loop ChPT plays an important role in understanding cur- 
rent LQCD results. However, numerical computations are not yet at a level of 
precision that 0(p 4 ) LECs can be reliably extracted from the lattice. The most 
extensive (small quark masses, large volumes) computations are done in the 
(partially) QA, to which ChPT can be adapted systematically. The fundamen- 
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Fig. 5. K + —+ tt + tt° decay amplitude as a lunction of pion mass, in units of the 
experimental value. Figure provided by JLQCD (Ishizuka et al. (1997)). Also shown 
are the data from Bernard and Soni (1989). Tree-level corrected. 



tal difference between QCD and its quenched relative shows up in full force at 
one-loop in ChPT: the nonanalytic terms in the QA are in general very different 
from those of the full theory. The safest way to look at these one-loop differences 
is to take them as an indication of the systematic error made by using the QA. 

We have seen that the quenched theory is afflicted with infrared divergences 
which show up as a diverging chiral limit, and a diverging infinite volume limit. 
This is all due to the special role of the r/ in the QA. It is therefore important to 
test the predictions of QChPT against lattice computations, as long as we will 
be using the (partially) QA. 
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